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Abstract 

We analyze the 1 + 1 dimensional Nambu-Jona-Lasinio model non- 
perturbatively. In addition to its simple ground state saddle points, the effective 
action of this model has a rich collection of non-trivial saddle points in which the 
composite fields a(x) = {jptp) and tt(x) = (^75^) form static space dependent 
configurations because of non-trivial dynamics. These configurations may be 
viewed as one dimensional chiral bags that trap the original fermions ( "quarks" ) 
into stable extended entities ("hadrons"). We provide explicit expressions for 
the profiles of these objects and calculate their masses. Our analysis of these 
saddle points is based on an explicit representation we find for the diagonal 
resolvent of the Dirac operator in a {a(x),7r(x)} background which produces a 
prescribed number of bound states. We analyse in detail the cases of a single as 
well as two bound states. We find that bags that trap N fermions are the most 
stable ones, because they release all the fermion rest mass as binding energy and 
become massless. Our explicit construction of the diagonal resolvent is based 
on elementary Sturm-Liouville theory and simple dimensional analysis and does 
not depend on the large N approximation. These facts make it, in our view, 
simpler and more direct than the calculations previously done by Shei, using 
the inverse scattering method following Dashen, Hasslacher, and Neveu. Our 
method of finding such non-trivial static configurations may be applied to other 
1 + 1 dimensional field theories. 
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1 Introduction 



Over the last thirty years or so, physicists have gradually learned about the behavior of 
quantum field theory in the non-perturbative regime. In 1 + 1 dimensional spacetime, 
some models are exactly soluble Another important approach involves the large 
N expansion f2|. In particular, in the mid-seventies Dashen, Hasslacher, and Neveu 
H used the inverse scattering method ||J] to determine the spectrum of the Gross- 
Neveu model ||. Recently, one of us developed an alternative method, based on the 
Gel'fand-Dikii equation ||, to study the same problem ]7] as well as other problems 
|| . As will be explained below, we feel that this method has certain advantages over 
the inverse scattering method. 

In this paper we study the 1 + 1 dimensional Nambu-Jona-Lasinio (NJL) Qmodel|J 
which is a renormalisable field theory defined by the action || 



A' 



S = (fx <Y1 ^aW^a + 



a=l 



N 



N 



\a=l I \a=l I 

describing N self interacting massless Dirac fermions ip a (a = I, . . . , N). This action 
is invariant under SU(N)f (g) U(l) ® U(1)a, namely, under 



and 

We rewrite ( |1 . 1|) as 
S = 



U ab ip b 



U e SU{N) 



f ' 



d 2 x {ip i$ — (a + i7T7 5 ) ip — — ^ (cr 2 + n 2 ) 



1 

27 2 



1.2) 



1.3) 



where a(x),7r(x) are the scalar and pseudoscalar auxiliary fields, respectively which 
are both of mass dimension 1. These fields are singlets under SU(N)f ® U(l), but 

1 This model is also dubbed in the literature as the "Chiral Gross- Neveu Model" as well as the 
"Multiflavor Thirring Model" . 

2 From this point to the end of this paper flavor indices are suppressed. Thus iip0ip should be 



N 



N 



understood as i^^^p a $ipa- Similarly ipTip stands for ip a Ttj} a , where T = 1. 



75- 



a=l 



transform as a vector under the axial transformation in (|1.2|), namely 



a + 275 7r — > e 2 * 75 ^(<t + 27571") 



Thus, the partition function associated with (|1.3|) is 



Z 



T>a Vtx Dip Dip exp i / d 2 x I ip 



(cr + iir^) 



if, 



2r 



7T 



:i-4) 



;i.5) 



Integrating over the grassmannian variables leads to Z — J T>oT>ix exp{iS e ff[a, tt]} 
where the bare effective action is 

1 



S eff [a, tt] 



2<T 



rf 2 x (a 2 + 



TT 



— iN Tr In 



— (cr + «VT7 5 ) 



1.6) 



and the trace is taken over both functional and Dirac indices. 

This theory has been studied in the limit N — ► oo with Ng 2 held fixed |J. In this 
limit ( |1.5|) is governed by saddle points of ( |1.6| ) and the small fluctuations around 
them. The most general saddle point condition reads 



SS P 



5a (x, t) 



5S 



8tt (x, t) 



a (x, t) 



71 (x, t) 
g 2 



+ iN tr 



N tr 



(x, t\ 



75 (x,t\- 



(CT + 27T7 5/ 



i0 - (a + Z7T7 5 ) 



\x, t) 



\x,t) 



0. 
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In particular, the non-perturbative vacuum of fll.l ) is governed by the simplest 
large N saddle points of the path integral associated with it, where the composite 
scalar operator ifjif) and the pseudoscalar operator itpj^tp develop space time indepen- 
dent expectation values. 

These saddle points are extrema of the effective potential V e ff associated with 
(|1.3|), namely, the value of —S e ff for space-time independent a, n configurations per 
unit time per unit length. The effective potential V e jf depends only on the combina- 
tion p 2 = cr 2 + 7r 2 as a result of chiral symmetry. V e ff has a minimum as a function 
of p at p = m ^ that is fixed by the (bare) gap equation 

d 2 k 1 



m + iNg 2 tr / 



(27T) 2 - 







m 
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which yields the dynamical mass 



m = Ae . (1.9) 



Here A is an ultraviolet cutoff. The mass m must be a renormalisation group invariant. 
Thus, the model is asymptotically free. We can get rid of the cutoff at the price of in- 
troducing an arbitrary renormalisation scale [i. The renormalised coupling gn (/i) and 

7T 1 J 

the cut-off dependent bare coupling are then related through A e N 3 2 w = fie N9 * M 
in a convention where Ng 2 R {m) = -. Trading the dimensionless coupling g\ for 
the dynamical mass scale m represents the well known phenomenon of dimensional 
transmutation. 

The vacuum manifold of (|1.3|) is therefore a circle p = m in the a, 7r plane, and 
the equivalent vacua are parametrised by the chiral angle 9 = arctan^. Therefore, 
small fluctuations of the Dirac fields around the vacuum manifold develop dynamical^, 
chiral mass mexp(i#7 5 ). 

Note in passing that the massless fluctuations of 9 along the vacuum manifold 



decouple from the spectrum ]T0[ so that the axial U(l) symmetry does not break 



dynamically in this two dimensional model JTTfl. 



Non-trivial excitations of the vacuum, on the other hand, are described semiclassi- 
cally by large N saddle points of the path integral over ( |1 . 1| ) at which a and 7r develop 
space-time dependent expectation values fl2|, [Tj[ . These expectation values are the 



space-time dependent solution of (|L7| ). Saddle points of this type are important also 
in discussing the large order behavior flU], |T5|j of the expansion of the path integral 



over (|1 . 1| ) . 

Shei [|TB| has studied the saddle points of the NJL model by applying the inverse 
scattering method following Dashen et al.||. These saddle points describe sectors of 
( |1 . 1| ) that include scattering states of the (dynamically massive) fermions in ( |l.ip , as 
well rich collection of bound states thereof. 

These bound states result from the strong infrared interactions, which polarise the 



3 Note that the axial U(l) symmetry in fll.2| ) protects the fermions from developing a mass term 
to any order in perturbation theory. 
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vacuum inhomogeneously, causing the composite scalar ipif) and pseudoscalar iijj^i/j 
fileds to form finite action space-time dependent condensates. These condensates are 
stable becuse of the binding energy released by the trapped fermions and therefore 
cannot form without such binding. This description agrees with the general physical 
picture drawn in ||17||. We may regard these condensates as one dimensional chiral bags 



18 , |19|| that trap the original fermions ( "quarks" ) into stable finite action extended 



entities ("hadrons"). 

In this paper we develop further the method of [[/], || , applying it to the N JL model 
(|1 . 1|) as an alternative to the inverse scattering investigations in [|l6j. We focus on 
static extended configurations providing explicit expressions for the profiles of these 
objects and calculate their masses. Our analysis of these static saddle points is based 
on an explicit representation we find for the diagonal resolvent of the Dirac operator 
in a <j(x), 7r(x) background which produces a prescribed number of bound states. This 
explicit construction of the diagonal resolvent can actually be carried out for finite 
N. It is based on elementary Sturm-Liouville theory as well as on simple dimensional 
analysis. All our manipulations involve the space dependent scalar and pseudoscalar 
condensates directly. In our view, these facts make the method presented here simpler 
than inverse scattering calculations previously employed in this problem because we 
do not need to work with the scattering data and the so called trace identities that 
relate them to the space dependent condensates. Our method of finding such non- 
trivial static configurations may be applied to other two dimensional field theories. 

It is worth mentioning at this point that the NJL model (|l.lj) is completely inte- 
grable for any number of flavors^] N. Its spectrum and completely factorised S matrix 
were determined in a series of papers |2IJ by a Bethe ansatz diagonalisation of the 
Hamiltonian for any number N of flavors. The large N spectrum obtained here as 
well as in |16| is consistent with the exact solution of |HJ. Note, however, that the 
large N analysis in this paper concerns only dynamics of the interactions between 



4 For N = 1 a simple Fierz transformation shows that ( |l . l| ) is simply the massless Thirring model, 
which is a conformal quantum field theory having no mass gap. A mass gap appears dynamically 
only for N > 2. 



4 



fermions and extended objects. We do not address issues like scattering of one ex- 
tended object on another, which is discussed in the exact analysis of |2(J. Consistency 
of our approximate large N results and the exact results of |2U| reassures us of the 
validity of our calculations. 

Rather than treating non-trivial excitations as abstract vectors in Hilbert space, 



which is inevitable in ||20|| , our analysis draws almost a "mechanical" picture of how 
"hadrons" arise in the NJL model. This description of "hadron" formation 
suit of inhomogeneous polarisations of the vacuum due to strong infrared interactions 
may have some restricted similarity to dynamics of QCD in the real world. Fur- 
thermore, our resolvent method is potentially applicable for non-integrable models in 
1 + 1 dimensions. In contrast, Bethe ansatz and factorisable S matrix techniques are 
limited in principle to 1 + 1 dimensions because of the Coleman-Mandula theorem |2l||, 
whereas large N saddle point techniques may provide powerful tools in analysis of 
more realistic higher dimensional field theories p2||. 

If we set n(x) in ( |1.3| ) to be identically zero, we recover the Gross-Neveu model, 
defined by 



a 2 



S GN = J d'x |V [iP ~ <t\ V - ^2 j ■ (1-10) 

In spite of their similarities, these two field theories are quite different, as is well- 
known from the field theoretic literature of the seventies. The crucial difference is 
that the Gross-Neveu model possesses a discrete symmetry, a — > —a, rather than the 
continuous symmetry ( |1.2| ) in the NJL model studied here. This discrete symmetry is 
dynamically broken by the non-perturbative vacuum, and thus there is a kink solution 
P3] , |3|, |7|, the so-called Callan-Coleman-Gross-Zee (CCGZ) kink cr(x) = mtanh(ma:), 
interpolating between ±m at x — ±oo respectively. Therefore, topology insures the 
stability of these kinks. 

In contrast, the NJL model, with its continuous symmetry, does not have a topo- 
logically stable soliton solution. The solitons arising in the NJL model and studied 
in this paper can only be stabilised by binding fermions. To stress this observation 
further, we note that the spectrum of the Dirac operator of the Gross-Neveu model 
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in the backgound of a CCGZ kink has a single bound state at zero energy, and there- 
fore no binding energy is released when they trap (any number of) fermions. The 
stability of the kinks in the Gross-Neveu model is guaranteed by topology already. 
In contrast, the stability of the extended objects studied here is not due to topology, 
but to dynamics. 

This paper is organised as follows: In Section 2 we prove that the static con- 
denstates cr(x) and tt(x) in (ft. .3 ) must be such that the resulting Dirac operator is 
reflectionless. Our proof of this strong restriction on the Dirac operator involves basic 
field theoretic arguments and has nothing to do with the large N approximation. We 
next show in Section 3 that if we fix in advance the number of bound states in the 
spectrum of the reflectionless Dirac operator, then simple dimensional analysis deter- 
mines the diagonal resolvent of this operator explicitly in terms of the background 
fields and their derivatives. We then construct the resolvent assuming the backround 
fields support a single bound state in Subsection 4.1. We are able to determine the 
profile of the background fields up to a finite number of parameters: the relative chiral 
rotation of the two vacua at the two ends of the one dimensional space and the bound 
state energies. In Subection 4.2 we provide partial analysis of the case of two bound 
states. We stress again that our construction of these background fields has nothing 
to do with the large N approximation. 

In order to determine these parameters we have to impose the saddle point condi- 
tion. We do so in Sections 5.1 amd 5.2. The relative chiral rotation of the asymptotic 
vacua is proportional to the number of fermions trapped in the bound states, in 
accordance with |24| . 



Some technical details are left to two appendices. In Appendix A we derive the 
spatial asymptotic behavior of the static Dirac operator Green's function. In order 
to make our paper self contained we derive the Gel'fand-Dikii equation in Appendix 
B. 
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2 Absence of Reflections in the Dirac Operator 
With Static Background Fields 



As was explained in the introduction, we are interested in static space dependent 
solutions of the extremum condition on S e g. To this end we need to invert the Dirac 
operator 



in a given background of static field configurations a(x) and n(x). In particular, 
we have to find the diagonal resolvent of (|2.1| ) in that background. The extremum 
condition on S e jf relates this resolvent, which in principle is a complicated and gen- 
erally unknown functional of <j(x), tt(x) and of their derivatives, to cr(x) and 7r(x) 
themselves. This complicated relation is the source of all difficuties that arise in any 
attempt to solve the model under consideration. It turns out, however, that basic 
field theoretic considerations, that are unrelated to the extremum condition, imply 
that fl2.1| ) must be reflectionless. This spectral property of Q2.1D sets rather pow- 
erful restrictions on the static background fields a(x) and n(x) which are allowed 
dynamically. In the next section we show how this special property of allows 
us to write explicit expressions for the resolvent in some restrictive cases, that are 
interesting enough from a physical point of view. 

Inverting ( |2.1| ) has nothing to do with the large N approximation, and conse- 
quently our results in this section are valid for any value of N. 

Here a(x) and tt(x) are our static background field configurations, for which we 
assume asymptotic behavior dictated by simple physical considerations. The overall 
energy deposited in any relevant static a, tt configuration must be finite. Therefore 
these fields must approach constant vacuum asymptotic values, while their derivatives 
vanish asymptotically. Then the axial U(l) symmetry implies that a 2 + ^ 2 x ~^ oc m 2 , 
where m is the dynamically generated mass, and therefore we arrive at the asymptotic 




(2.1) 
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boundary conditions for a and 7r, 



a — > mcos#± 

a:— >±oo 



a' ► 

>±oo 



msin#-i 



7T > 7/H51±l(7-|- 

x— >±oo 



7T' ► 

>±oo 



(2.2) 



where 6± are the asymptotic chiral alignment angles. Only the difference 9 + — 6_ is 
meaningful, of course, and henceforth we use the axial symmetry to set 6L = 0, such 
that er(— oo) = m and 7r(— oo) = 0. We also omit the subscript from 9 + and denote 
it simply by 9 from now on. It is in the background of such fields that we wish to 
invert Q2.1|) . 

In this paper we use the Majorana representation 7 = 02 , 7 1 = «cr 3 and 7 5 = 



-77 



G\ for 7 matrices. In this representation ( |2.1| ) becomes 



/ -a 



D 



G —ILO — 17T 



\ iuj — m d x — o 
Inverting ( |2.3| ) is achieved by solving 

/ — d x — a(x) — iuj — iix{x) \ ( a(x,y) b(x,y) 



(2.3) 



\ iuj — m{x) 



d x - a(x) 



-ilS(x-y) (2.4) 



V c(x,y) d(x,y) 

for the Green's function of ( |2.3| ) in a given background a(x),n(x). By dimensional 
analysis, we see that the quantities a, b, c and d are dimensionles. 

The diagonal elements a(x,y), d(x,y) in (27 3) may be expressed in term of the 
off-diagonal elements as 



a(x,y) 



i [d x - a(x))c(x,y) 



d(x,y) = i[d * + aix) }\ {x > y) (2.5) 



uj — ir(x) " lu + n(x 

which in turn satisfy the second order partial differential equations 

b(x,y) 



-d T . 



d x b(x, y) 




to + ir(x) 


+ 


d x c(x, y) 


+ 


UJ — 7r(x) 



a{x f + n{x f _ a > {x) _ w 2 + cr(x)n'(x) 

uj + 7r(x) 



a(x) 2 + tt(x) 2 + a'(x) - lu 2 + 



a{x)ir [x) 

UJ — 7l(x) 



UJ + 7t(x) 

c(x,y) 

UJ — 7r(x) 



5(x - y) 

-5(x - y) . 
(2.6) 



Thus, b(x, y) and —c(x, y) are simply the Green's functions of the corresponding 
second order Sturm-Liouville operators in Q2.6|) , 

6(x — y) b 2 (x)b 1 (y) + 9 (y - x) b 2 (y)b 1 (x) 



b(x,y) 



W h 



, v {x - y) c 2 {x)c x (y) + 6 (y - x) c 2 (y)c x {x) 
c{x } y) = — . (2.7) 



Here b\(x) and b 2 (x) are the Jost functions of the first equation in (|2.6| ) and 



w = b 2 (x)b\(x) -b!(x)b' 2 (x) 
b uj + tt(x) 

is their Wronskian. The latter is independent of x, since b x and b 2 share a common 
value of the spectral parameter uj 2 . Similarly, Ci,c 2 are the Jost functions of the 
second equattion in (|2.6|) and W c is their Wronskian. We leave the precise definition 



of these Jost functions in terms of their spatial asymptotic behavior to Appendix 
A, where we also derive the spatial asymptotic behavior of the static Dirac operator 
Green's function. Substituting ( |2.7| ) into ( |2.5| ) we obtain the appropriate expressions 
for a(x,y) and d(x,y), which we do not write explicitly. 5 

We define the diagonal resolvent (x |i_D _1 |x ) symmetrically as 



(x | — iD \x) 



A{x) B(x) 
C{x) D(x) 



- lim 

2 e-*0+ 



/ a(x, y) + a(y, x) b(x,y) + b(y, x) 



(2.9) 

\ c(x,y) + c(y,x) d(x,y) + d(y,x) 
Here A(x) through D(x) stand for the entries of the diagonal resolvent, which follow- 
ing ( p.5|) and (|2.7p have the compact representational 

bi(x)b 2 (x) _ i [d x + 2a(x)\ B (x) 



B(x) 



W h 



UJ + 7r(x) 



5 It is useful however to note, that despite the d x operation in (2.5), neither a{x,y) nor d(x,y) 
contain pieces proportional to S(x — y) . Such pieces cancel one another due to the symmetry of 
(2.7) under x <-> y . 

A, B, C and D are obviously functions of uj as well. For notational simplicity we suppress their 
explicit uj dependence. 
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C(x) 



ci(x)c 2 (x) 



A{x) 



i [d x - 2a(x)} C (x) 



(2.10) 



2 lo — ir(x) 

A simplifying observation is that the two linear operators on the left hand side of 
the equations ( |2.6| ) transform one into the other under a simultaneous sign flip[] of 
a(x) and tt(x). Therefore c(a,n) = —b(—a, —n) , and in particular 

C(a,n) = -B(-a,-n), (2.11) 

and thus all four entries of the diagonal resolvent ( |2.9| ) may be expressed in terms of 
B(x). 

The spatial asymptotic behavior of (|2.9|) is derived in Appendix A and given by 



(|A.6|) . A more compact form of that result is 

1 + R(k)e 2ik ^ r 



(x | — iD 1 \x) 



x— >±oo 



2k 



i7 5 7r(x) — a(x) — -y°u; 



+ 



R(k)e 



2ik \ x\ 



7 1 sgnx 



(2.12) 



where k = yju 2 — m 2 and R(k) is the reflection coefficient of the first equation in 
Note that for u 2 > m 2 , i.e., in the continuum part of the spectrum of (|2.3|), the 



piece of the resolvent Q2.12|) that is proportional to R(k) oscillates persistently as a 
function of x. This observation has a far reaching result that we now derive. Consider 
the expectation value of fermionic vector current operator j M in the static cr(x), tt(x) 
background^ 



(a(x), Tr{x)\f\a(x), ir(x)) 



du 
— 1 
2tt 



-tr 



7" 



A(x) B(x) 
C(x) D(x) 



(2.13) 



7 This is merely a reflection of the fact that coupling the fermions to ^75 does not respect charge 
conjugation invariance. 

8 In the following it is enough to discuss only the vector current, because the axial current j£ = 
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Therefore, we find from fl2.12|) that the asymptotic behavior of the current matrix 
elements is 



(a(x), 7r{x)\j \a{x), tt(x)) — > 

x— >±oo 

and 

(a(x), 7r(a;)|jV(x), tt(x)) _> - / ^i?(A;)e 2ifc ^ (2.14) 

x— >±oo J Z7T 

where we used the fact that / Y~i:f{k) = because k (u) is an even function of uj. 

Thus, an arbitrary static background a(x) , tt(x) induces fermion currents that do 
not deacy as x — > ±oo, unless R(k) = 0. Clearly, we cannot have such currents in 
our static problem and we conclude that as far as the field theory ( |1.3| ) is concerned, 



the fields a(x) , ir(x) must be such that the Sturm-Liouville operators in ( |2.6| ) and 
therefore the Dirac operator (|2.3|) are reflectionless. 

The absence of reflections emerges here from basic principles of field theory, and 
not merely as a large N saddle point condition, as in ||, [l(J. Indeed, reflection- 
lessness of ( |2.3| ) must hold whatever the value of N is. Therefore, the fact that 
reflect ionlessness of ( [2.3|) appeared in [||, |16| as a saddle point condition in the inverse 



scattering formalism simply indicates consistency of the large iV approximation in 
analysing space dependent condensations a(x) ,tt(x). The absence of reflections also 
restores asymptotic translational invariance. What we mean by this statement is that 
if R(k) = then ( |2.12j ) is simply the result of inverting ( |2.3|) in Fourier space with 
constant asymptotic background ( |2.2|) , namely, 



(x 



iD^lx) 



( imcosO uj + msin9 



2 - ,,c 



2 



(2.15) 



y — uj + msin# imcos9 
which therefore yields the asymptotic behavior of ( |2.9|) for properly chosen chiral 
alignment angles. Note that in the absence of reflections, ( |2.12| ) attains its asymptotic 



value ( p,15| ) by simply following the asymptotic behavior of cr(x) and of n(x), which 
are the exclusive sources of any asymptotic x dependence of the resolvent. This 
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expression (|2.15| ) has cuts in the complex u plane stemming from scattering states 



of fermions of mass m. These cuts must obviously persist in A, B, C and D away 
from the asymptotic region, and we make use of this fact in the next section. We 
used the asymptotic matrix elements (|2.13|) of the vector current operator in the 
background of static cr(x), tc(x) to establish the absence of reflections in the static 
Dirac operator. We can now make use of this result to examine its general dynamical 
implications on matrix elements of other interesting operators, namely, the scalar ipip 
and pseudoscalar tpi^if) density operators. Their matrix elements in the background 
of <t(x), ii(x) are 



and 



(<j(x), ir(x)\il>ip\a(x), tt(x)} = N I — tr (x\iD 1 \x) 



(<y(x), 7r(x)\il)iy5ip\a(x), tt(x)) = N I — tr ry 5 (xliD 1 \x) . (2.16) 



Therefore, from (|2.12|) their asymptotic behavior is simply 

duj l + R(k)e 2ik ^ 



(cr(x), ir(x)\i/jilj\a(x), vr(x)) — > Na(x) 



x^±oo J 2tt k 

and 

r du 1 + i? (k) e 2ik \ x \ 
(a(x), n(x)\^ l5 ^\a(x), Nn(x)) — vr(x) / ^ V • (2.17) 

x^±oo J Z7T K 

Clearly, in the absence of reflections, the asymptotic x dependence of these matrix 
elements follows the profiles of <j(x) and tt(x), respectively. Otherwise, if R(k) ^ 0, 
these matrix elements will have further powerlike decay in x superimposed on these 
profiles, which is not related directly to the typical length scales appearing in <y(x) 
and in tt(x). We close this section by investigating implications of ( |2.17| ) for extremal 



background configurations. For such configurations the matrix element of the scalar 
density is equal to a(x)/g 2 and that of the pseudoscalar density is equal to n(x)/g 2 . 
Such background fields must obviously correspond to a reflectionless Dirac operator, 
but let us for the moment entertain ourselves with the assumption that R(k) in (|2.17|) 



is arbitrary and see how the absence of reflections appears as a saddle point condition. 
Thus, for extremal configurations, as x —>■ ±oo, <j{x) cancels off both sides of the first 

12 



equation in ( |2.17|) and tt(x) cancels off both sides of the other equation. This leaves 
us with a common dispersion integral 



du 1 + R (k) e 



2ik \x\ 



Ng 2 J 2tt k 

It turns out (see (|5.21| below) that the integral over the first x independent term on the 
right hand side cancels precisely the constant term on the left hand side. This is simply 
a reformulation of (|1.8| ) in Minkowski space. Therefore, the remaining x dependent 
integral must vanish for any (large) It follows then, that R(k) must vanish. 

Thus absence of reflections appears here as a saddle point requirement, in a rather 
simple elegant manner, without ever invoking the inverse scattering transform. The 
whole purpose of this section is to prove that one cannot consider static reflectionful 
backgrounds to begin with, and thus the emergence of absence of reflections as a 
saddle point condition is simply a successful consistency check for the validity of the 
large N approximation applied to space dependent condensates. 



13 



3 The Diagonal Resolvent for a Fixed Number of 
Bound States 



The requirement that the static Dirac operator ( |2.3| ) be reflectionless is by itself quite 
restrictive, since most a(x),7[(x) configurations will not lead to a reflectionless static 
Dirac operator. Construction of explicit expressions for the resolvent in terms of 
a(x) ,tt(x) and their derivatives is a formidable task even under such severe restric- 
tions on these fields. We now show how to accomplish such a construction at the 
price of posing further restrictions on a(x) and ir(x) in function space. However, even 
under these further restrictions the results we obtain are still quite interesting from 
a physical point of view. 

In the following we concentrate on the B(x) component of (p.9| ). The other entries 
in ( |2.9| ) may be deduced from B(x) trough ( |2.10|) and ( [2.11| ). 

Our starting point here is the observation that one can derive from the represen- 
tation of B(x) in ( |2.10| ) a functional identity in the form of a differential equation 
relating B(x) to cr(x) and tt(x) without ever knowing the explicit form of the Jost 
functions bi (x) and 6 2 (x). We leave the details of derivation to Appendix B, where 
we show that the identity mentioned above is 



d r 



u> + ir(x) 



d. r 



' d x B(x) ' 

UJ + 7l(x) 



UJ + 7r(x) 

2B(x) 
[uj + tt(x)] 



\d x 


B(x) 




) 




UJ + 7r(x) 





a{xf + n(x) 2 - a'(x) - uj 2 + < x >) x ] 

UJ + Tl{X) 



- 3 

2 Ux 



cr(x) 2 + 7r(x) 2 — cr'(x) + 



<t(x)7t'(x) 



= 



(3.1) 



uj + ir(x) 

with a similar expression for C(x) in which a — > —a it — > —it that we do not write 
down explicitly. 

Here we denote derivatives with respect to x either by primes or by partial deriva- 
tives. This equation is a linear form of what is referred to in the mathematical 
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literature as the "Gel'fand-Dikii" identity 0. This identity merely reflects the fact 
that B{x) is the diagonal resolvent of the Strum-Liouville operator discussed above 
and sets no restrictions on cr(x) and ir(x). 

If we were able to solve (|3.1|) for B(x) in a closed form for any static configuration 
of <j(x),tt(x), we would then be able to express (x|i-D _1 |x) in terms of the latter 
fields and their derivatives, and therefore to integrate (|1.7|) back to find an expression 
for the effective action (|1.6| ) explicitly in terms of a(x) and n(x). Invoking at that 
point Lorentz invariance of ([L.6|) we would then actually be able to write down the full 
effective action for space-time dependent a and 7r. Note moreover that in principle 
such a procedure would yield an exact expression for the effective action, regardless 
of what N is. 

Unfortunately, deriving such an expression for B(x) is a difficult task, and thus we 
set ourselves a simpler goal in this paper, by determining the desired expression for 
B(x) with cr(x), tc(x) restricted to a specific sectors in the space of all possible static 
configurations. To specify these sectors consider the Dirac equation associated with 
(|2.1|) , Dip = 0. For a given configuration of a(x), ir(x) (such that D is reflectionless), 
this equation has n bound states at energies Ux,...,uj n as well as scattering states. 
A given sector is then defined by specifying the number of bound states the Dirac 
equation has.| 

As we saw above, B(x) must have a cut in the u plane with branch points at 
uj = ±to. If in addition to scattering states a(x),ir(x) support n bound states at 
energies u\, ...0J n (which must all lie in the real interval — m < uj < m)Q then the 
corresponding B must contain a simple pole for each of these bound states. Therefore, 



9 The effect of scattering states on B(x) is rigidly fixed by spatial asymptotics, as ( 2.15 ) indicates, 
so only bound states are used to specify such a sector. 

10 The Gross-Neveu model[|[ ||, |?j is a theory of Majorana (real) fermions. Therefore its spectrum 
is invariant under charge conjugation, i.e., it is symmetric under u> — ► —u>. Thus in that case the 
bound states are paired symmetrically around u> — and B(x) is really a function of oj 2 . The chiral 
NJL model on the other hand is a theory of Dirac (complex) fermions, charge conjugation symmetry 
of the spectrum is broken by the 7r field and bound states are not paired. 
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B(x) must contain the purely u dependent factor 

1 



(3.2) 



Vm 2 -u 2 Uk=i(^ - u k) 
of mass dimension — n — 1. Any other singularity B(x) may have in the complex uo 
plane cannot be directly related to the spectrum of the 

Dirac operator, and therefore must involve x dependence as well. Based on our 
discussion in Appendix A, the only possible combination that mixes these variables 
is exp(iv^ 2 — nT' 2 x). But such a combination is ruled out as we elaborated in the 
previous section, by the requirement that the Dirac operator be reflectionless. The 
factor ( |3.2j ) then exhausts all allowed singularities of B(x) in the complex u plane. 
Recall further that B(x) is a dimensionless quantity, and thus the negative dimension 
of the uo dependent factor ( |3.2[ ) must be balanced by a polynomial of degree n + 1 in 
uj (with x dependent coefficients) of mass dimension n + 1, namely^ 

^ _ B n+1 {x)u n+1 + .... + B^uo + Bq{x) 
Vm 2 -w 2 IlLi( w -^) 

The mass dimension of Bh(x) (k — 0, n + 1) is n + 1 — k. 

The main point here is that simple dimensional analysis in conjunction with the 
prescribed analytic properties of B(x) fix its u dependence completely, up to n + 1 
unknown bound state energies, and n + 2 unknown functions of a(x), ir(x) and their 
derivatives. These functions are by no means arbitrary. They have to be such that 
( |3.3| ) and the resulting expression for C (x) are indeed the resolvents of the appropriate 
Sturm-Liouville operators. These expressions for B(x) and C(x) must be therefore 
subjected to the Gel'fand-Dikii identities ( |3.1| ) and the corresponding identity for 
C{x). 

Substituting B(x) into fl3.1|) we obtain an equation of the form 



Q$ 5 (cu,x)/[u; + ir(x)} 4 = 0, (3.4) 



11 One may argue that Eq. ( |3.3[ ) should be further multiplied by a dimensionless bounded function 
/(— )• However such a function must be entire, otherwise it will changed the prescribed singularity 
properties of B(x), but the only bounded entire functions are constant. 
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where Q„ +5 (u, x) is a polynomial of degree n + 5 in u with x dependent coefficients 
that are linear combinations of the functions Bk(x) and their first three derivatives. 

Note that because of the linearity and homogeneity of ( |3.1|) , the purely u> depen- 
dent denominator of ( |3.3|) with its explicit dependence on the bound state energies 
drops out from ( |3.4j ). This is actually the main advantage[^]of working with the linear 



form of the Gel'fand-Dikii identity rather than with its non-linear form ( B.10 )[f7|. 

Setting to zero each of the x dependent coefficients in Q n + 5 we obtain an overde- 
termined system of n + 6 linear differential equations in the n + 2 functions Bk(x). 
Using n + 2 of the equations we fix all the functions Bk{x) in terms of cr(x), ir(x) and 
their derivatives, up to n + 2 integration constants bf.. These integration constants 
are completely determined once we enforce on the resulting expression for B(x) the 
asymptotic behavior fl2.2p and ( |2.15| ). The integration constants bk turn out to be 
polynomials in m 2 and the bound state energies u^. 

At this stage we are left, independently of n, with four non-linear differential 
equations in cr(x) and 7r(rr)0. A similar analysis applies for C(x), leading to an 
equation of the form (w,x) = 0, where following Q2.11|) (u),a(x),ir(x)) = 

~~Qn+5 \ u i ~ cr ( x ): ~ TC ( X )) ■ Setting the first n + 2 coefficients in Q„ +5 (uj,x) to zero 



we verify that C(x) is related to B(x) as in ( 2.11 ), but that the remaining four 

equations for <j(x) and tt(x) are different from their counterparts associated with 

B(x) as the explicit relation between Q^+ 5 and Q^+b suggests. We are thus left with 

an overdetermined set of eight non-linear differential equations for the two functions 

a{x) and n(x). Observing that Q„+ 5 ± Qn+b is °dd (even) in a and it, we note that 

these eight equations are equivalent to breaking each of the four remaining equations 

associated with B(x) into a part even in a and tt and a part odd in a and it and 

setting each of these parts to zero separately. 

Mathematical consistency of our analysis requires that the six most complicated 

12 The non-linear version of the Gel'fand-Dikii identity (B.8) (or ( B.10| )) contains further informa- 
tion about the normalisation of B(x) 7 but the latter may be readily determined from the asymptotic 
behavior ( J2.15Q of B(x). 

13 Coefficients of the various terms in these equations are also polynomials in m 2 and the bound 
state energies u>k- 
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equations of the total eight be redundant relative to the remaining two equations, 
because we have only two unknown functions, cr(x) and tt(x). This requirement must 
be fulfilled, because otherwise we are compelled to deduce that there can be no <r(x) 
and tt(x) configurations for which the Dirac equation Dip = has precisely n bound 
states, with n = 0, 1, 2, • • •, which is presumably an erroneous conclusion. 

Therefore cr(x) and tt(x) are uniquely determined from the two independent equa- 
tions given the asymptotic boundary condittions (|2.2f) they satisfy. This leaves only 
the bound state energies undetermined, but the latter cannot be determined by the 
resolvent identity, which does not really care what their values are. These energies are 
determined by imposing the saddle point conditions fll.7|), i.e., by dynamical aspects 
of the model under investigation. 

In the preceding paragraphs we laid down the mathematical aspects of our anal- 
ysis. We now add to these a symmetry argument which will simplify our solution 
of the differential equations for <r(x) and ir(x) a great deal. The two non- redundant 
coupled differential equations for cr(x) and tt(x) allow us to eliminate one of these 
functions in terms of the other. We choose to eliminate^] ir(x) in terms of a(x), 

ir a (x) = G a [a a (x)} (3.5) 

where a is a global chiral alignment angle. This relation is clearly covariant under 
axial rotations a — > a + Aa, because a(x) and tt(x) transform as the two components 
of a vector under U(1)a as ( |1.4|) shows. We expect (|3.5|) to be a linear relation. 
Imposing the boundary conditions ( |2.2|) we have 

Q 

ir(x) = —[a(x) — m] c °t- • (3.6) 

In this way we reduce the problem into finding the single function a(x). The 
condition ( |3.6|) is an external supplement to the coupled differential equations for 
a(x) and tt(x) stemming from the Gel'fand-Dikii equation. We thus have to make 
sure that the resulting solution for a(x) and ( |3.6j ) are indeed solutions of these coupled 
differential equations. 

14 We prefer to eliminate n(x) in terms of a(x) because the latter never vanishes identically. 
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We now provide the details of such calculations in the case of a single bound state, 
as well as partial results concerning two bound states. 



19 



4 Extended Object Profiles 
4.1 A single bound state 

In this case ( |3.3| ) becomes 

= B 2 (x)uj 2 + Bi(x)u + B Q (x) ^ 
(uj — uji) y/m 2 — uj 2 

where the single bound state energy is uj%. Then setting to zero the coefficients of uj 6 
through cj 4 in the degree six polynomial (|3.4j ), we find 

B 2 (x) = b 2 , 



Bi(x) = b 2 7r(x) + b\ and 



.80(0;) = biir(x) + -j- \o 2 (x) + n 2 (x) - a\x)} + b 



(4.2) 



where b 2 , b\ and b Q are integration constants. We then impose the asymptotic bound- 
ary conditions 

1 



B(x) 



— uj' 



(uj + msm9±) 



to fix the latter, 



bi 



' 2 



m 
T 



and therefore 



B(x) 



uj + ir(x) ^ a 2 (x) + 7r 2 (x) — a (x) — m 2 



2y/m 2 — uj 



4 (uj — UJi) \J Tfl 2 — UJ 2 

The relation ( |2. 1 1|) then immediately leads to 

UJ — 7r(x) 



C(x) 



cr 2 (x) + ir 2 (x) + cr'(x) — m 2 



2^/m 2 — uj 2 



(4.3) 



(4.4) 



(4.5) 



4 (u — UJ\) \/m 2 — uj 2 
Having the coefficients of uj 6 through uj 4 in the degree six polynomial ( |3.4|) set to 
zero, we are left with a cubic polynomial 



43. 



(jn 2 — T[ 2 (x) — o 2 (x) s j n(x) — ujio'(x) H — ix (x) 



+ 



UJi (n 2 (x) + a 2 (x)^j - a(x)ir'(x) + a' (x)n(x) } uj 3 H = 



(4.6) 
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where the (■ • •) stand for lower powers of id. The cubic ( |4.6|) has to be set to zero 
identically, producing eight coupled differential equations in o{x) and u{x) as we 
discussed above. The simplest equation of these is obtained by setting to zero the 
part of the a; 3 coefficient in ( f4.6|) that is even in a and it, namely, 



Or. 



u>i (tt 2 (x) + <J 2 (x)^j — <j(x)it (x) + a (x)n(x) 







which we immediately integrate into 



tt 2 (x) + cr 2 {x) — m 2 = <7(x)tt (x) — a (x)tt{x) 



(4.7) 



Here we have used the boundary conditions ( |2.2| ) to determine the integration 
constant. The next simplest equation is obtained by setting to zero the part of the 
u 3 coefficient in ( [4.6| j that is odd in a and it, and so on. 

Following our general discussion we solve the system of coupled equations (|3.6|) 
and ([4.71) , which leads to 



d_ 

dx 

Solving (fOP we find 



a[x) — m 



2lu, 



2uj\ tan 7 



ct(x) — m msm.6 



m sin9 tan 7 



<7[X) 



7c(x) 



m 



1 + exp 2o; 1 tan| • (x — x ) 
msin8 



(4.9) 



1 + exp 2cj 1 tan2 • (x — x ) 

where we have chosen the integration constant (parametrised by xo) such that <j(x) 
and 7r(x) would be free of poles. Note that the boundary conditions at x — ► +oo 
require 



9 

u i tan- < . 
2 



(4.10) 



Substituting the expressions (O) into (O) one finds that the resulting B(x) is 
indeed a solution of the corresponding Gel'fand-Dikii equation ( p.l|) , verifying the 
consistency of our solution. Our results (|4.9| ) for a(x) and 7r(x) agree with those of 
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|]IE|| . They have the profile of an extended object, a lump or a chiral "bag", of size 
of the order cot|/k;i centred around an arbitrary point x$. Note that the profiles in 
f) satisfy 



p 2 (x) 



a 2 (x) +7T 2 (x) 



m 



m 2 sin 2 (#/2)sech 2 



e , 

c^tan- • [x — Xq) 



(4.11) 



Thus, as expected by construction, this configuration interpolates between two differ- 
ent vacua at x = ±00. As x increases from —00, the vacuum configuration becomes 
distorted. The distortion reaches its maximum at the location of the "bag", where 
m 2 — p 2 (xq) = m 2 sin (6/2) and then relaxes back into the other vacuum state at 
x = 00. The arbitrariness of Xq is, of course, a manifestation of translational invari- 
ance. 
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4.2 Two Bound States 

In this case (13.3) becomes 



_ B 3 (x)u 3 + B 2 (x)u 2 + B 1 (x)u + Bpjx) 
(c<j — ooi) [uj — w 2 J V m — ^ 

where the bound state energies are u)\ and co>2- (Obviously, -B(x) in this subsection 
should not be confused with its counterpart in the previous subsection.) In this case 
the polynomial (|3.4j ) is of degree seven in uj. Following the procedure outlined in 
Section 3 we find after imposing the boundary conditions fl2.2| ) that 

v. LU + Tc(x) 

B(x) = — v ' + 
2\ / m 2 — uj 2 

[a 2 (x) + 7r 2 (x) — a (x) - m 2 } [n (x) + u - uj\ - uj 2 ] 
4 (lj — U)i) {uj — 0J2) \Jm 2 — uj 2 

tx" (x) - 20-0)71-' 0) 



8 (uj — u)\) (uj — UJ2) \Jm 2 — io 2 
Then, from ( 2.11|) we find that 

C(x) - w - * {X) 



(4.13) 



2y/m 2 — uj 2 

a 2 (x) + tc 2 (x) + a (x) — m 2 } [—-k(x) + u — uo\ — uj 2 ] 
4 (uj — ujx) (uj — UJ2) \Jm 2 — ui 2 

tx" (x) + 2a(xW (x) 



(4.14) 



8 (a; — tux) (uj — uj 2 ) y/m 2 — uj 2 
Note that if we set uii +uj 2 = and tc(x) = the resolvents ( |4.13|) and ( |4.14| ), and 



therefore the whole spectrum, become invariant under uj — > —uj, and we obtain the 
equation appropriate to the Gross-Neveu model. 

Setting to zero the coefficients of uj 7 through uj 4 in the degree seven polynomial 
( |3.4j ), we are left with a cubic polynomial in u which we do not write down explic- 
itly. This polynomial must vanish identically, producing eight coupled differential 
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equations in a(x) and tt(x) as we discussed above. The simplest equation of these is 
obtained by setting to zero the part of the uo 3 coefficient in that polynomial which 
unlike the previous case, is now odd in a and tc 

2d x { 4 (wi + uj 2 ) [tt 2 (x) + a 2 {x) - m 2 } tt(x) + 6 [vr 2 (x) + a 2 (x)} a'{x) 

+2 (2ujiuj2 — m 2 ) cr'(x) — 2 (u)\ + uj 2 ) ^ (x) — er (x) } = 0. (4.15) 

As in the previous case, this is a complete derivative which we readily integrate into 

4 (wi + w 2 ) [tt 2 (x) + a 2 (x) - m 2 ] vr(x) + 6 [tt 2 (x) + cr 2 (x)] cr'(x) 

+2 (2cj!CJ 2 - m 2 ) a'(x) - 2 (wi + w 2 ) 7r"(a?) - </"(a;) = . (4.16) 



Here we have used the boundary conditions ( |2.2j ) to determine the integration con- 
stant. Note that ( 4.16|) is of third order in derivatives and cubic, whereas its single 



bound state counterpart (|4.7| ) is only first order in derivatives and quadratic. 

The next simplest equation is obtained by setting to zero the part of the a; 3 
coefficient in the cubic polynomial that is evem in a and tt, 

2 d x { 2 (m 2 - 2lo!U 2 ) [tt 2 {x) + a 2 {x)] - 3 [vr 2 (x) + <x 2 (x)] 2 + 4(wi + u 2 ) 
■[a(x)7r'(x)-a'(x)7r(x)]}+4[7r(x)7r'"(x) + a(x)a'"(x)\ =0. (4.17) 

and so on. 

Following our general discussion we have to solve solve the system of coupled 
equations ( |3.6| ) and ( |4. 16| ) which is equivalent to 



2 A (wi + uj 2 ) [4m y 2 + 2(1 + \ 2 )y 3 - y"} 

+d x { 4 (m 2 + W1W2) y + 6 m y 2 + 2 (A 2 + 1) y 3 - y"} = (4.18) 

where A = — cot(6 l /2) and y(x) = cr(x) — m. We have not succeeded in solving this 
non-linear ordinary differential equation in closed form. 



24 



5 The Saddle Point Conditions 



Derivation of the explicit expressions of ct(x) and tt(x) does not involve the saddle 
point equations (|1.7|). Rather, it tells us independently of the large N approximation 
that a(x) and ir(x) must have the form given in ( |4.9| ) in order for the associated 
Dirac operator to be reflectionless and to have a single bound state at a prescribed 
energy uo\ in addition to scattering states. Thus, for the solution (|4.9|) we have 
yet to determine the values of U\ and 9 allowed by the saddle point condition. More 
generality, our discussion in Section III will lead us to the cr(x) and ir(x) configurations 
which correspond to reflectionless Dirac operators with a prescribed number of bound 
states at some prescribed energies u>i,u>2, ■•■ in addition to scattering states. As 
emphasized earlier, this result is independent of the large A" limit. The allowed 
values of 9, a>i,a>2, • • • must then be determined by the saddle point condition (|1.7|) . 
It is this dynamical feature that we can analyse only in the large A" limit. 

For static background fields the general saddle point condition flL7 ) assume the 
simpler form 

du 



a(x)+ Ng 2 f ^[A(x) + D(x)] = 
J Ztt 



ir(x) + iNg 2 [ ^ [B(x) + C{x)\ = 0. (5.19) 
J 2n 

In Subsection 5.1 we impose this condition on the explicit single bound state back- 
ground we found in the previous section and calculate the mass of such "bags" . The 
two bound states case is discussed in Subsection 5.2. 

5.1 A single bound state 

Substituting ( |4.4j ), ( |4.5| ) into the saddle point equations (|5.19|) we obtainF] 
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In the following formula wc omit explicit x dependence of the fields. The number 7r also appears 
in the formula, but only in the combination ^ . Therefore there is no danger of confusing the field 
tt(x) and the number w. 
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a f duo 

+ i 



Ng 2 
+2 



2tx A\/m 2 — uo 2 (uj — u)\) (uo 2 — n 2 ) 



{Aouo* + 2 (tt' - 2uj x a) uo 2 



11 1\ i L n 

a [a —7i — m J — Co>i7r — -a 



tt f duo 2ttu — 2oj-i 7r — a' 

+ i 1 



Ng 



-A 



2n 2\Jm 2 — u 2 {uj — uj\) 



uo - 2n 2 (tt' - 2ui<t) 



0. 



(5.20) 



These equations are dispersion relations among the various x dependent parts 
and u)\. Clearly, both dispersion integrals in ( 5.20| ) are logarithmically divergent in 
A, but subtracting each of the integrals once we can get rid of these divergences. 
The required subtractions are actually already built in in ( p.20| ). In order to see this 
consider the (bare) gap equation ( |1.8| ) in Minkowski space. This equation is equivalent 
to the (logarithmically divergent) dispersion relation^ 



1 



2n \Jm 2 — uj 2 + ie 



(5.21) 



If we now replace each of the coefficients in ( |5.20| ) by the integral on the right 
hand side of ( [5 -21 ) 0, the divergent parts of each pair of integrals cancel and the 
equations ( |5.20| ) become 



duo tt' + 



^F(a, tt) 



2tx \lw? — uo 2 (uo — UJ\) 



duo 



a 



2ix yjm 2 —uj 2 (uj — uj\) 



where 



a,7T 



a(a 2 + 7r 2 - m 2 ) - uo x tx' a" 



and C is the contour in the complex uo plane depicted in Fig.(l) 



(5.22) 



(5.23) 



3 To see this equivalence simply perform (contour) integration over spatial momentum first. 
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m 



CO 



777 



Fig. 1: The contour C in the complex u> plane in Eq. ( 5.22) ) . The continuum states appear as 
the two cuts along the real axis with branch points at ±m (wiggly lines), and the bound state is the 
pole at u>. 



The expression ( |5.23p is the residue of the x dependent poles at lu = ±tt(x) 
in the first equation in ( [5.221 ). The quantisation condition ( p.22[ ) on uj\ cannot be x 
dependent. Therefore fl5.23|) must vanish as a consistency requirement. Substituting^ 
(p|) in ( F23|) we find that 

. , mtan 2 t 
F(x) = 



sech 



9 

u>i xtan- 
2 



cos#e~ Wl xta,n 2 -)_ e Wl xtan 2 



2 2 2^ 

to-, — m cos - 
. 1 2 



Thus, F((T, it) vanishes for the configurations ( |4.9|) provided 



1 2 

m cos 



6. 



(5.24) 



which sets an interesting relation between the bound state energy and the chiral 
alignment angle of the vacuum at x — > +oo. This relation actually leaves 9 the only 
free parameter in the problem with respect to which we have to extremise the action. 
The condition ([4.1 0Q then picks out one branch of ( |5.24| ). 

We still have to determine U\ and 9 separately. Now the saddle point condition 
simply boils down to the single equation 

du 1 



2n ^/m? — uj 2 {uj — uj\) 



0. 



(5.25) 



r Here we have set xq = for simplicity. 
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The contour integral in (|5.25| ) is most conveniently calculated by deforming the 
contour C into the contour C shown in Fig. (2). The "hairpin" wing of C picks up 
the contribution of the filled Fermi sea, and the little circle arount the simple pole at 
uj = uoi is the contribution of fermions populating the bound state of the "bag" . 



Itn CO A 



- m 




CO 



Rett) 



m 



Fig. 2: The deformed integration contour C leading to Eq. (5.26) 



Assuming that the "bag" traps n fermions in that state, and recalling that each 
state in the Fermi sea continuum accomodates N fermions, we see that Q5.25| ) becomes 

iN 



mr — ui 



7T 



m + u}\ n 

arctanw h — — 1 

V m — uj\ N 







(5.26) 



The solution of this equation yields the quantisation condition 



mcos 



mi 



(5.27) 



in agreement with [I? 

It follows from (pj), ( p^|) and (gjg) that 

27m 



9 



N 



(5.2? 



that is, the relative chiral rotation of the vacua at ±oo is proportional to the number 
of fermions trapped in the "bag" . Our result ( p.28[ ) is consistent with the well known 
result that the fermion number current in a soliton background can be determined 
in some cases by topological considerations p4|- Note from (|5.28[) , that in the large 
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N limit, 9 (and therefore u{) take on non-trivial values only when the number of the 
trapped fermions scales finite fracion of N . 

As we already mentioned in the introduction, stability of "bags" formed in the NJL 
model are not stable because of topology. They are stabilised by releasing binding 
energy of the fermions trapped in them. To see this more explicitly, we calculate now 
the mass of the "bag" corresponding to ( |4.9|) , (|5.27|) and (|5.28|) . The effective action 
(|1.6|) for background fields ( [4.9|) is an ordinary function of the chiral angle 9 0. Let 
us denote this action per unit time by S(9)/T. Then, from (|2.9|) and ( |5.21|) we find 



1 OS 
T~d9 



dxduj 
2vr 



+i 



ia[x) 



ir(x) 



(A + D) 



da(x) 
89 



or 



[B + C) 



dir(x) 
09 



(5.29) 



Then as in our analysis of the saddle point condition (which is simply the condition 
0) we use (fOj), (jL5|), (plO|) and the fact that ( |5T23D vanishes to find 



for — 



Then, 



i_ds_ 

T~d9 



dx dui d x a(x)dgir(x) — d x 7i(x)dga(x) 



2n 2 (uj — uji) y/m 2 — u 2 

leads to the factorised expression 



]_dS_ 
T~d9 



4sin z 



dx (a — m)a' / — 



c 2ir (a; — oji) \Jm 2 — uj 2 



(5.30) 



where C is the contour in Fig. (2). The space integral is immediate and is essentially 
fixed by the boundary conditions (|2.2|) . The spectral integral is given by the left hand 



side of ( 5.26| ), but with a generic loi given by to\ = mcos(|). Here we have chosen 



the particular branch of (|5.24j ) that contains all the extremal values of uoi. Putting 
everything together we finally arrive at 



1 dS 

Wr~d9 



m 
~2 



n_ JT 
N + 2vr 



9 



sin- 



(5.31) 



18 Recall that io\ is a function of and not a free parameter. 
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The zeros of ( |5.31| ) are simply the zeros of ( |5.26| ), as these two equations are one and 



the same extremum condition. Integrating ( |5.31| ) with respect to 9 we finally find 

_1 -S{9) = (± + f] cos 9 - - hm 9 ~. (5.32) 



NTm y ' \N 2n J 2 vr 2 
Note that Q5.32j ) is not manifestly periodic in 9 because the Pauli exclusion principle 
limits 9 to be between and 27r. 

The mass of a "bag" containing n fermions in a single bound state is given by 
—S/T evaluated at the appropriate chiral angle ( |5.28| ). We thus find that this mass 
is simply 

Nm ixn 

M n = sin — 5.33 

7T N 



in accordance with 16, 20 1. These "bags" are stable because 



. n{n 1 +n 2 ) irn, . im 2 , K Q/n 

sin < sm + sin 5.34 

N N N K ' 

for ni,n 2 less than N, such that a "bag" with n\ + n 2 fermions cannot decay into two 
"bags" each containing a lower number of fermions. 

Entrapment of a small number of fermions cannot distort the homogeneous vac- 
uum considerably, so we expect that M n will be roughly the mass of n free massive 
fermions for n << N. As a matter of fact we used this expectation to fix the inte- 
gration constant in ( |5.33| ). For n « N we have M n ~ nm[l — | (jjf^J +■■•], so the 
binding energy released 

B -~-(iv) + "' (5 - 35) 

is indeed very small. However, as the number of fermions trapped in the "bag" 
approaches N, M n vanishes and the fermions release practically all their rest mass 
Nm as binding energy, to achieve maximum stability ]T7|]. In a weakly coupled field 
theory containing solitons, the mass of these extended objects is a measure of -\, 
the inverse square of the coupling constant. Here we have 4^ = N. It is amusing to 
speculate that these maximally stable massless solitons may teach is something about 
the strong coupling regime of the NJL model. 
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Note from ( |5.28| ), that the soliton twists all the way around as the number of 
fermions approaches N. In this case uj± — > —m, and the pole the resolvent has at 
uj = uj\ pinches the branchpoint uj = —m at the edge of the filled Dirac sea. One may 
wonder whether this enhanced singularity is a mathematical artifact, as the bound 
state simply tries to plunge into the filled Dirac sea. But this is clearly not the case. 
Indeed, uj\ is occupied by N fermions (in a flavor singlet). Their common spinor wave 
function must still be part of the discrete spectrum of the Dirac operator, because 
the highest lying state of the sea at uj = — m is already occupied by a flavor singlet 
made of N fermions, sharing a continuum spinor wave function, and therefore Pauli's 
exclusion principle protects the bound state from "dissolving" into the sea. 

5.2 Two bound states 



We concluded Subsection 4.2 short of an explicit solution of (|4.18| ), namely, short of 
an explicit expression for the two bound state background fields a(x) and n(x). In 
the following we make the eminently reasonable assumption that such a background 
exists, and pursue our analysis of its saddle point condition as far as we can without 
having its explicit form in hand. 

As in the previous subsection, we substitute ( f4.13| ) and ( (4.14| ) into the saddle point 
equations ( |5.19|) . We then make use of ( 5.21|) to eliminate the ultraviolate logarithmic 
divergences and to write the saddle point conditions as 



C" 

duj 



2^ 



ia(x) 



\/m 2 — 



uj" 



(A + D) 



K + uj 



U) 2 —1T 2 (x) 



2n 2\Jm 2 — m 2 (to — u>i) (to — uj 2 ) 



rmand 
du 
2^ 

C" 

duj 



\/m 2 — 



(B + C) 



M + uja' (x) 



I 2ir 2\/m 2 — m 2 {uj — u)\) {uj — co> 2 ) 



(5.36) 
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where C" is a contour similar to the contour C in Fig. (2) that encircles the additional 
pole at u 2 as well, and K(x), L(x) and M(x) are given by 



K(a,<ir) 



-a(a 2 + 7T 2 - m 2 ) + {u x + cu 2 )n' + \-o" 



(id! + L0 2 ) 



( 1 i 2 2> 



-0" ' 



V 2 + 7r 2 — m 2 



+ Co>iCg> 2 + cr 2 7t' 



and 
M(a,n) 



7T 

+ r 



-7r(a 2 + 7t 2 — m 2 ) — (cui + 002)^ + - 71 "" 



(5.37) 



Note that K (cr, 7r) differs from —F(a, tt) in (|5.23|) only by the additional term 
co> 2 vr'. The expression L(cr, tt) is the residue of the x dependent poles at uj = ±tt(x) in 
the first equation in ( |5.36| ). The quantisation conditions ( |5.36| ) on u>i and u; 2 cannot 
be x dependent. Therefore L(a, tt) must vanish as a consistency requirement. As we 
do not have the explicit expressions of a(x) and 7r(x), we assume from now on that 
L indeed vanishes. This is the only extra assumption we make. Then, assuming that 
the "bag" traps ni fermions in ui and n 2 fermions in u 2 , ( [5.36 ) boils down to the 
simple conditions 



J(wi) = I(u 2 ) 
where I(oo) is given in ( [5.261 ). Therefore, 



m cos 



which are identical in form to single bound state energy levels. From the general 



m cos 



fn 2 7r 



(5.38) 



(5.39) 



considerations of pj| we expect that the chiral angle 9 will be proportional to the 
total number of fermions trapped by the "bag" , so (|5.28 ) must read now 



e 



2-k (m + n 2 ) 
N 



(5.40) 
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The soliton mass is a function of m and of the chiral angle 9. Assuming this function 
is the same as in the previous case we therefore conjecture that the mass of the two 
bounnd state "bag" is simply 

M ni , m = — sm^ '- . (5.41) 

As far as mass is concerned, such a "bag" cannot be distinguished from a single bound 
state "bag" containing the same total number ri\ + n 2 of trapped fermions. If our 
conjecture is true, then such "bags" are stable against decaying into several "bags" 
with lower numbers of fermions as ( |5.34j ) shows. 
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A Appendix A 



In this Appendix we provide precise definitions of the Jost functions bi through c 2 in 
terms of their spatial asymptotic behavior and derive the spatial asymptotic behavior 
of the static Dirac operator Green's function. 



We concentrate for the moment on the first equation in ( |2.6|) . The boundary 
conditions (|2.2|) lead to the following simple spatial asymptotic behavior 



[-d 2 x + m 2 - u 2 } b(x) = 

of the homogeneous part of that equation. Thus, solutions of that homogeneous 
equation assume the generic asymptotic form 

M+e ikx + N+e -ikx ^ x _^ +oc 

b(x,u)~{ (A.l) 



M _ e ikx + N _ e -zkx X ^ -OO 



where 



k(u) = +Vlo 2 - m 2 . (A.2) 
On the real u axis k(u) is real for \u\ > m, which corresponds to scattering states of 



2j|) . Bound states of ( |2.3| ) reside in the domain \u\ < m, where k(u) = +i\/m 2 — u 2 = 
+in{uj) is purely imaginary and lies in the upper half plane. 

The Jost functions bi and b 2 alluded to in Section 2 form a particular pair of 
linearly independnt solution of the homogeneous equation mentioned above, specified 
by their asymptotic behavior. Let the asymptotic amplitudes of b r (x), (r = 1, 2) in 



([A. 1|) be M r ±,N r ±. The asymptotic form ( A.l ) of b\(x) has by definition Mi_ = 



and that of b 2 {x) has = 0. One may summarise our definitions of b\ and b 2 , by 
saying that b\ corresponds to a one dimensional scattering situation where the source 
is at +oo emitting waves to the left (the term N 1+ e~ thx ) and that b 2 corresponds to 
a one dimensional scattering situation where the source is at — oo emitting waves to 
the right (the term M 2 _e lkx ). Note also that outside the continuum, b\(x) decays to 
the left while b 2 (x) decays to the right. With these definitions the Wronskian (|2.8|) 
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becomes 



W b (+oo) = -2ik- 



M 2+ N : 



UJ + TC ( + OOJ 



-2ik- 



M 2 _N 1 _ 

UJ + 7T ( — OO) 



W h {-oo) . (A.3) 



Therefore, it follows from (|2.10 ) and from ( 2.11| ) that the entries A,B,C, and D in 
(|2.9|) have the asymptotic form 



A{x) 
B(x) 
C(x) 

D(x) 
clS X — > ±00, where 



■^-{[a(x) -iksgnx}R{k)e 2iklxl + a(x)} 



UJ + 7r(x) " 

-2ik 

UJ — 7t(x) 

2ik 



1 + R(k) e 



2ik\x\ 



1 + R(k) e 



2ik\x\ 



^-\\a(x) + iksgnx}R(k)e 2lklxl + <r(x)\ (A.4) 
2k 1 > 



R{k) 



M 1+ N 2 



(A.5) 



N 1+ M 2 - 

is the reflection coefficient of the Sturm-Liouville operator in the first equation in 



( 2.6|) . The diagonal resolvent of the Dirac operator is therefore 



ia(x) uj + ty(x) 



-UJ + 7r(x) ia(x) 



+ 



iR(k)e 



2ik 



| x .| / io (x) + k sgncc uj + ti(x) 



2k 



-uj + 7r(x) zct(x) — sgnx 



(A.6) 
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B Appendix B 



Consider the Sturm-Liouville problem 

— p(x)if) (x) + \V(x) — Ep(x)]if)(x) = 



-oo < x < oo . 



(B.l) 



We assume that the "metric" p(x) does not vanish anywhere and that the weight 
function p{x) is positive everywhere. E is a complex number, called the spectral 
parameter. 

As in our discussion in the main text and in the previous appendix, let ipi{x) be 
the Jost function which decays as x — ► — oo for values of E below the continuum 
threshold. Similarly, let ip2(x) be the Jost function which decays as x — > +00. Then, 
the Green's function of the operator in (|B.1|) is 



G{x,y) 



6{x-y) ^ 2 (g)^i(y) + Q (y - x) j) 2 {y)^i{x) 
W 



(B.2) 



where 



W = p(x) ■0 2 (^)'0i( :E ) ~ 4 J l{ x )' l l>2( X ) 



(B.3) 



is the (x independent) Wronskian of these two functions. Note that (|B.2|) decays (at a 
rate dictated by the Jost functions) as either one of its argument diverges in absolute 
value, holding the other one finite, as long as E does not hit one of the eigenvalues 
of the Sturm-Liouville operator. 

As in the main text the diagonal resolvent R(x) = G(x,x) is defined as 



R(x) 



- lim [G{x, x + e) + G(x + e, x)} = r v Z K ' 



We then use O and (PH) to show that 

ipi pR' + 1 tp' 2 pR' — 1 

ih ~ 2pR ' fa ~ 2pR 



(B.4) 



(B.5) 



Finally, using ( |B.1| ) and ( |B.5| ) we find 
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(pR*)' = 2(V - Ep)R + 



(pR' 



2pR 



(B.6) 



and 



\p{pR')'\ = [2p(V - Ep)R]' + 2p(V - Ep)BL . (B.7) 

Note that the non-linearity of (P-6| ) in R has disappeared after one more differen- 
tiation with respect to x. 

Multiplying (|B.6|) through by 2R we find 



- 2pR(pR')' + {pR') 2 + ApR 2 (V — pE) = 1 (B.8) 

is the linear form of the Gel'fand- 



which is the Gel'fand-Dikii equation [y]. Eq. 
Dikii equation we use in the text (Eq. ([O] .) 

The quantities corresponding to the discussion of the Dirac operator in the text 

are 



p(x) = p(x) 
and 

V(x) ' 



to + ir(x) 



E = 00 



a(x) 2 + n(x) 2 - a'(x) + 



to + n(x) 

The Gel'fand-Dikii equation ( |B.8| ) then reads 



(t(x)tt'(x) 

LO + n(x) 



(B.9) 



B{x) n 
2 -t-t9 x 

CO + 7T(x) 



d x B{x) 

LO + k{x) 



+ 



d x B(x) 

CO + Tl(x) 



n 2 



+ 



2B{x) 

00 + 7r(x) 



1 2 



a(x) 2 + n(x) 2 - a\x) - oo 2 + lXj 



i 



(B.10) 
> 0. 



oo + 7r(a;) 

Strictly speaking, Sturm-Liouville theory requires that p(x) = p(x) 
Our solution for 7r(x) turns out to be bounded, so all formulae are valid a posteriori 
for large positive oo. Such a restriction on oo, though mathematically required, is 
unphysical. 
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Note however, that because of the relation (|2.11 ), we may view C(x) as a contin- 
uation of B{x) to large negative uj. 

An important application of the Gel'fand-Dikii identities ( |B.7| ), ( |B.8| ) is that they 
generate an asymptotic expansion of R in negative odd powers of y/~E. The explicit 
uj (and therefore E) dependence of our specific p(x),p(x) and V(x) complicates this 
expansion. 
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